Linear magnetoresistance in commercial n-type silicon due to inhomogeneous doping 
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Free electron theory tells us that resistiv- 
ity is independent of magnetic field^. In fact, 
most observations match the semiclassical predic- 
tion of a magnetoresistance that is quadratic at 
low fields before saturating^. However, a non- 
saturating linear magnetoresistance has been ob- 
served in exotic semiconductors such as silver 
chalcogenides^'^, lightly-doped InSb^, N-doped 
InAs^, MnAs-GaAs composites!, PrFeAsO§, and 
epitaxial graphene^. Here we report the obser- 
vation of a large linear magnetoresistance in the 
ohmic regime in commonplace commercial n-type 
silicon wafer. It is well-described by a classical 
model of spatially fluctuating donor densities^", 
and may be amplified by altering the aspect ratio 
of the sample to enhance current-jetting: increas- 
ing the width tenfold increased the magnetoresis- 
tance at 8 T from 445 % to 4707 % at 35 K. This 
physical picture may well offer insights into the 
large magnetoresistances recently observed in n- 
type-*^-*^ and p-type^^ Si in the non-ohmic regime. 

The conventional theory of magnetoresistance (MR) 
in metals and semiconductors relies upon a distribution 
of scattering times amongst the conducting carriers that 
cannot be compensated by a unique Hall field^. For a 
material with a closed free electron Fermi surface (FS) 
and a principal charge carrier this leads to positive MR 
that is quadratic in weak magnetic fields, B, and satu- 
rates in strong magnetic fields according to: 
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where p is the resistivity. In the denominator of equa- 
tion[T]it can be shown^^that CB^ = (Zmf/?'c)^ = (^ct)^, 
where ujc is the cyclotron frequency, r the mean free time 
between scattering events, Imf the mean free path and 
the cyclotron radius. In the strong field limit {uJcT ^ 1), 
electrons complete several orbits before scattering. In the 
weak field limit (wcT ^1), only a fraction of these orbits 
are completed, and so Ap/p ^ AB^, depending predom- 
inantly upon the anisotropy in relaxation times and cy- 
clotron masses of the carriersi^. Other galvanomagnetic 
effects are present in magnetically ordered materials^. 

In the past decade some materials have emerged that 
break these previously established trends. Narrow gap 
semiconductors such as iron diantinomide^^, the silver 
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chalcogenides^, indium antimonido^. and nitrogen-doped 
indium arsenide^ have shown unexpectedly high MR over 
large temperature ranges. In each case, the MR was 
non-saturating, approximately linear at high fields, and 
was enhanced by strong disorder arising from impurity 
substitutioni^iii and spatial inhomogeneity of the sam- 
ple stoichiometry^ii^. In the last few years numerical 
modelling of large scale inhomogeneous inclusions using 
classical models has been used to describe this linear 
MRi^^. 

Here we have measured the MR in n-type commercial 
Si wafer, which we attribute to a more subtle mecha- 
nism due to the weak disorder from spatially fluctuat- 
ing conductivity that arises statistically from the ran- 
dom donor distributioriF^. Si:P wafers with 3.6 ± 0.3 
51cm room temperature resistivity and dopant density 
(1.4 ± 0.1) X 10^^ cm""^ were used. The net carrier den- 
sity, (n), and an absolute value of the carrier mobility, 
/i, were obtained as a function of temperature, T, using 
the van der Pauw method on square wafers (see supple- 
mentary figure 1). The MR was measured using four in- 
line contacts on a sample with the following dimensions: 
width, w = 1.8 ± 0.1 mm; length, I = 5.5 ± 0.5 mm; 
thickness, t — 530 ± 50 /zm; and voltage probe separa- 
tion, s — 2.0 ± 0.4 mm. Indium ohmic contacts provided 
a four wire current- voltage I-V characteristic that was 
ohmic at low bias throughout all temperatures. Mea- 
surements were performed in three orientations of field 
with respect to current density j: transverse (j _L B_) 
with the magnetic field out of the sample plane; perpen- 
dicular (j ± B) with the magnetic field in plane; and 
longitudinal (j \\ B_). 

The transverse MR in the ohmic regime in figure [TJi 
was approximately linear up to high magnetic fields and 
showed no indication of saturation in 8 T fields at tem- 
peratures from 30 - 200 K, and is much larger than at 
room temperature where the system barely enters the 
high-field regime^. The high field gradient of this plot 
is known as the Kohler slope (KS). The data were re- 
duced to a Kohler plot of Ap/ p versus WcT, shown in[Tl3. 
In weak and strong fields the MR scaled onto a unique 
curve according Kohler's law; Ap/p = /(uJct), where / 
is a scaling function, using the relationship lOcT = fj.B. 
Longitudinal measurements did not obey Kohler's law 
but possessed a variable KS as a function of temperature 
as shown in the supplementary figures 2 and 3. In weak 
fields (shown in the inset of figurefl})) the departure from 
quadratic behaviour (shown as dashed line) occurred for 
less than half a completed cyclotron orbit, while linear 
MR was satisfied in strong fields for WcT > 2. 

Linear MR can arise in polycrystallinc material with 
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Figure 1: Transverse MR at: 30 (■), 40 (•), 50 (a), 60 (t), 
80 {■<), 120 (►) and 200 K (♦) in transverse magnetic fields, 
a) MR increased monotonically with decreasing temperature. 
In tlie strong magnetic field condition MR was non-saturating 
and approximately linear, b) MR at all temperatures reduced 
down to one unique curve provided nB — ljcT conforming to 
Kohler's rule. Tlie data is magnified in the inset to demon- 
strate the weak to strong field crossover at ljcT ~ 1 and uni- 
versal scaling of the weak field MR. A quadratic law is plotted 
at weak fields, locT <C 1 (dashed line) 



an open FS resulting from the averaging of saturating 
MR from closed orbits and non-saturating MR from open 
orbits2^. Such effects can also be caused by quantum 
linear magnetoresistancc (QLMR)^i2^. The former can 
be disregarded for the single crystals used here and the 
closed silicon conduction band energy surfaces. The 
latter requires a very small n and effective mass, m* 
which generally limits it to metals and semi-metals. As 
m* ~ 0.4me in Si:P— , a field of 10 T would not produce 
QLMR above 30 K— since the thermal energy would ex- 
ceed the separation of the Landau levels, ruling out this 
mechanism here. 

Aside from these possible causes, Pippard proposed 
four distinct mechanisms that could lead to strong field, 
non saturating linear MR in crystals^: small angle scat- 
tering, magnetic breakdown, voids, and inhomogeneities. 

The influence of small angle scattering is limited 
at higher temperatures, where electrons are scattered 
through larger angles. In Si:P at 100 K, where a large 
positive MR is observed, the most effective phonons 
would scatter at ^ 30° suggesting that small angle scat- 
tering is not applicable^. Equally, magnetic breakdown 
from tunnelling between orbits in momentum-space is un- 
likely as the FS of the conduction band in Si lies far from 
the Brillouin zone boundary. 

MR arising from spherical voids with a volume fraction 
F scales as OASGFuIcT in a transverse field and 0.637FuJcT 
in a longitudinal field^. As the magnitude of the trans- 
verse MR shown in figure[TJi is much greater than the lon- 
gitudinal MR (sec supplementary figure 2) the observed 
MR does not agree with this hypothesis. In figure [Ud, 
with the aspect ratio l/w ~ 3, the KS was 0.30 ± 0.01, 



implying a void volume fraction F = 0.62 ± 0.02 in our 
single crystal Si wafer, which is utterly implausible. 

Thus, linear MR from inhomogeneities remains as the 
most probable cause of enhanced MR in strong fields. 
The limit of strong disorder has been used to explain the 
linear MR in the silver chalcogenides2£. However, sili- 
con growth techniques aim to produce uniformly doped 
wafers and thus the disorder is likely to be weak. The 
suspicion that fluctuation in dopant density of pulled 
crystals may influence the conduction of semiconduct- 
ing single crystals was the original motivation for Her- 
ring's paper in 1960^^^ although there was little direct 
experimental evidence available at that time. The spa- 
tial variations in Hall coefficient that these defects give 
rise to were predicted to prevent saturation of the trans- 
verse MR whilst having no effect on the saturating longi- 
tudinal MR. For instance, in the non-degenerate weakly 
doped sample discussed here, the average impurity sep- 
aration is ~ 90 nm. The mean free path at 200 K, 
^MF = {fJ./e){3'm*kBTy^^, where e is the electron charge 
and fee is Boltzmann's constant, would be ~ 92 nm^l. 
As the electrons will scatter on a Icngthscale compara- 
ble to the impurity separation they can be susceptible 
to fluctuations in the concentration of ionised donors, 
A^i, on lengthscales greater than their mean free path. 
Herring suggested that the fractional variance (FV) of 
the number of ionised donors {{Ni ~ (Ai))^)/(A^i)^ was 
the relevant measure of the fluctuation magnitude. If 
this value was signiflcant over lengthscales greater than 
a suitably deflncd effective mean free path, ^off, and the 
Debye screening length. Ad, then the inhomogeneity that 
results from the distribution of random impurities would 
create uncompensated Hall fields that contribute signifi- 
cantly to a linear MRf^. 

The fractional variance should hence be limited on a 
lengthscale Xd = ^efcRT/ {{n)e'^), where e is the Si:P 
permittivity, in lightly doped samples, and by ^cff = 
(2/(a;cT)^)~^/'^ZMF in highly doped samples. Isotropic, 
saturating MR such as in equation [1] is assumed in the 
case of homogeneous doping. When donor density fluc- 
tuations are introduced, the effective MR for transverse 
isotropic fluctuations in can be written asi^: 

Apcff Ap B A{N,-{N,)f) 
Pes ^ P nep{Or (A^i)2 ' ^ ' 

where C has a numerical value that depends upon the 
conductivity anisotropy. In the isotropic case relevant 
here, C = tt/4. 

To estimate the FV of a random distribution of im- 
purities over the relevant lengthscales, we simulated a 
cubic crystal consisting of 10^^ atoms, substituting ran- 
dom host atoms for ionised impurities. The FV of the 
impurities probed over different lengthscales is shown in 
figure [2^. The average density of ionised impurities at 
50 and 200 K was taken from supplementary figure 1 as- 
suming that Ni n in ohmic conduction. Included in 
figure [2^. arc the lower limits on the FV at each tempera- 
ture given by the lengthscales Ad and l^ff, calculated for 
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Figure 2: Interpretation of the observed MR within the Her- 
ring model, (a) Simulation of the impurity concentration fluc- 
tuations of a cubic crystal of 10^^ atoms on different length- 
scales. The vertical lines mark the Debye screening length, 
Ad, and effective mean free path, Zeff, at 50 K and 200 K. 
These parameters provide an upper limit for the FV to which 
the electron will be sensitive. For this impurity concentra- 
tion the FV is limited by Zcff at both temperatures but still 
remains large, b) Fits (lines) to the transverse MR (symbols) 
at 50 and 200 K according to the simplified model of Herring 
in equation (2] 



= 2. These limitations suggest that the FV should 
be no more than 0.27 at 50 K and 0.77 at 200 K. 

To see if these values of the FV in impurity density 
are sufficient to produce the observed linear MR, fits to 
equation [5] were performed using the measured ^, /o(0), 
and (n) values with only the FV and A as fitting param- 
eters. The fits are shown in figure [^b, showing that a 
FV of 0.34 is required at 50 K and 0.21 at 200 K. At 50 
K, it is expected that the value for for which these 
conditions are satisfied should be 2; from our fits, the 
value for lo^t ~ 2.25 is in reasonable agreement. Overall, 
the classical Herring model provides a very satisfactory 
explanation for our observations. 

This non-saturating linear MR of the Si:P can be en- 
hanced by classically understood methods, by altering 
the aspect ratio of the wafers, which causes current jet- 
ting. When attempting to discern the physical MR of 
semiconductors, such geometrical enhancement is to be 
avoided, as above. However, it can be used to enhance 
the response of devices in magnetic fields. Solin et al. 
have shown that the geometrical enhancement in high 
mobility semiconductors in certain geometries such as 
Corbino disks can be used for practical applications such 
as read heads for high-density recording^^*^. 

We prepared nine further samples with increasing 
width, w, with indium contacts positioned at a separation 
I ~ 10 mm between current leads and s ~ 3 mm between 
voltage probes. Photographs of these samples, labelled 
A-H, are shown in figure[3]with the dimensions, measured 
using vernier calipers, given in the caption. Perpendic- 
ular and transverse MR at 35 K of samples A to H are 
shown in figures[2^ and[3]3 respectively. In the perpcndic- 
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Figure 3: Geometrical enhancement MR, measured in per- 
pendicular (a) and transverse (b) field orientations at 35 K 
shown mirrored to one another using the right hand axis to 
encourage comparison. Wafers A-H are shown with respective 
values of a = Ijw: 10 ± 1 (A,B), 4.5 ± 0.6 (B,«), 3.7 ± 0.5 
(C,A), 2.6±0.3 (D,T), 2.3±0.2 (E,^), 1.9±0.3 (F,^), 1.6±0.2 
(G,*), 1.1 ±0.2 (H,*). 



ular orientation, where the sample projection in the field 
direction does not change, the variation in MR between 
the samples is very small. In the transverse field shown 
in figure [31d the projection of the sample in the magnetic 
field direction varied significantly. In this orientation, in- 
creasing the width of the samples dramatically enhanced 
the MR of the films. Sample A, the narrowest sample, 
reached a maximum Ap/p = 445 % at 8 T. The widest 
sample, H, had a maximum Ap/p = 4707 %, over an 
order of magnitude larger. 

Geometrical enhancement of transverse MR has been 
calculated for rectangular plates of homogeneous conduc- 
tors with electrodes spanning the sample ends22i2^:our 
geometry differs slightly in that we have point-like con- 
tacts, but the same general conclusions hold. For the 
narrowest sample A, the net current flow would be ap- 
proximately homogeneous providing a good indication of 
the physical MR intrinsic to the material. In sample H, 
the widest sample, current flow would be inhomogeneous 
with significant current deflection perpendicular to the 
magnetic field near to the electrodes, increasing the volt- 
age detected there and hence the measured resistance. 

Hence, we conclude that the large non-saturating MR 
we observe here relies upon the distribution of ionised 
donors being sufficiently disordered that current distor- 
tions arise from the local fluctuating Hall fields. The 
donor distribution in this study is random, defined by a 
three-dimensional Poisson distribution, but with modern 
lithographic and ion implantation techniques, the local 
density could be tailored to maximise this effect. Alter- 
natively, for a thin sample, lithography techniques could 
produce an array of top gates that could tune the local 
density of ionised donors to provide strong disorder by 
locally modifying the Hall constant, providing a MR re- 
sponse that can be tailored with an applied electric field. 

If the carrier concentration is decreased, as is the case 
for very lightly doped silicon, the linear MR from the FV 
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in equation [2] is limited by the Debye electron screening 
length. This limitation from screening could perhaps be 
overcome by the large electric fields resulting from the 
concentration of conduction electrons. As Ad cx 
in the non-quasineutral Mott-Gurney regime, Ad would 
be reduced without changing the FV of the ionised impu- 
rities. As the FV would remain large, with the reduced 
limitation from Ad, a large FV may produce a large MR 
even in intrinsic samples. This could help to explain the 
large MR in high electric fields reported for intrinsic sili- 
con (i-Si)iii2ai^. For instance, with typically (iVj) - 10^^ 
cm3 at 300 K, {{N; - (iVi))2)/(7Vi)2 - 0.01. By increas- 
ing the source-drain current to inject more carriers such 
that (jiinj) ^ 10{n) this value would increase to ~ 0.5. 
It is possible that the significant enhancement that could 
arise in both Mott-Gurncy regime and during electrical 
breakdown could increase the MR to produce the large 
measured valuefM^. 



Methods 



In all wafers the (001) crystal vector was normal to the surface. 
High angle XRD verified that the wafers were cleaved along (110) 
orientations such that the current, j, was parallel to (110) in all 
measurements. 



Contacts were made to the wafer surface by mechanical cleavage 
through the native oxide with a diamond scribe followed by impres- 
sion of indium. The contacts were then annealed at 350 °C for 10 
minutes. When the contacts were made, aluminium wire was at- 
tached using a wedge bonder from the sample holder to the indium. 
A 3 X 3 mm'^ square wafer was used for van der Pauw measure- 
ments using magnetic fields of ±1.5 T to ascertain Hall voltages. 
All magnetoresistance measurements were performed in a gas flow 
cryostat using a four probe in-line DC technique at low bias in the 
ohmic regime, with the current bias adjusted to maintain roughly 
5 mV between the voltage probes. 
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